Draft version September 24, 2012 

Preprint typeset using I^TeX style emulateapj v. 5/2/11 



A THEORETICAL COLOR- VELOCITY CORRELATION FOR SUPERNOVAE ASSOCIATED WITH 

GAMMA-RAY BURSTS 

Sharon Rapoport^'^, Stuart A. Sim^'^, Keiichi Maeda^, Masaomi Tanaka^ ,Markus Kromer^, Brian P. Schmidt^'^, 

Ken'ichi Nomoto^ 

Draft version September 24, 2012 

ABSTRACT 

We carry out the first multi-dimensional radiative transfer calculations to simultaneously compute 
synthetic spectra and light curves for models of supernovae driven by fast bipolar outflows. These 
allow us to make self-consistent predictions for the orientation dependence of both colour evolution 
and spectral features. We compare models with different degrees of asphericity and metallicity and 
find significant observable consequences of both. In aspherical models, we find spectral and light 
curve features that vary systematically with observer orientation. In particular, we find that the early 
phase light curves are brighter and bluer when viewed close to the polar axis but that the peak flux 
is highest for equatorial (off-axis) inclinations. Spectral line features also depend systematically on 
observer orientation, including the velocity of the Si ll 6355A line. Consequently, our models predict a 
correlation between line velocity and color that could assist the identification of supernovae associated 
with off-axis jet-driven explosions. The amplitude and range of this correlation depends on the degree 
of asphericity, the metallicity and the epoch of observation but we find that it is always present and 
acts in the same direction. 

Subject headings: Gamma-ray burst :general, supernovaigeneral, radiative transfer 



1. INTRODUCTION 



Gamma-ray bursts (GRBs) provide a power ful tool to 
study galaxy formation at high redshifts (e.g. "Conselice] 
|et al 2005) and provide a unique opportuni ty to ob serve 
extreme physical proce sses (e.g. [Bucciantini et al.||2009 ; 
Woosley fc Heger||2006^ . Long duration GRBs (LGRB 



have been shown to b e connected with broad-lined Type 
Ic supernovae (SNe Ic lGalama et a l.|1998| Iwamoto et al.' 
■1998; Kulkarn iend1 ll998; Bloom erar| P99 ; Meszaros 
2006j [Lindner et al. 2010), implying that they are linked 
to the collapse of a massive star, see Woos ley fc Bloom] 



( [20061 ) for 

a review. However, the exact nature of this 



connection is not yet understood, and is especially poorly 
constrained in the very bright GRBs used to study the 
high-redshift universe. 
The collap sar model ([Woosley 1993 [MacFadyen fc 



Woosley|l999| ) predicts that when the iron core of a mas- 



sive rotatmg star collapses into a black hole (BH), mate- 
rial accretes onto the BH via a rotating disk. A narrow, 
highly relativistic jet is launched at the black hole - ac- 
cretion disk boundary which gives rise to the observed 
gamma-ray burst. In addition, a sub-relativistic wind 
is driven off the disk, which drives the explosion of the 
star. Nucleosynthesis in the explosion pr oduces ^^Ni, the 



decay of which powers the SN emission ( Iwamoto et al. 
[T998 ) . Angular momentum is believed to be a crucial lac- 
tor in determining the mass that collapses into the BH. 
A more massive BH implies lower density in the accre- 
tion disk, which might disfavour production of a GRB 
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(Woosley & Zhang 2007|). As higher metallicity leads 
to greater mass loss via stellar winds, which causes the 
star to lose angular momentum, GRB progenitors are ex- 
pected to have low metallicity (| Heger et al.||2003 



Yoon 



Langerl2005l). In the close binary progenitor sce nario, 



the m etallicity dependence is much weaker (e.g. Fryer 
(eFaL 2007) 

Observationally, this collapsar model implies that 
emission produced in the highly collimated jet is strongly 
affected by relativistic beaming such that it is only de- 
tectable when the Earth is closely aligned with the jet 
direction. Therefore, detection of a GRB immediately 
suggests that the event is observed pole-on (along the 
rotation axis). In contrast, the associated SN could be 
observed from all orientations, although its properties 
can be expected to depend on viewing-angle. 

Studying the SN associated with LGRBs is one means 
to better understand the explosion mechanism. Analysis 
of SN 1998bw, the first SN believed to host a LGRB, 
supports the high- mass collapsar model ([iwamoto et al. 



1998[ [Woosley et al.||1999| [MacFadyen <^ Woosley ||1999p 
and gives an estimated main sequence progenitor mass of 
~ 30 — 4OM0. However, this object's 7-ray luminosities 
were orders of magnitude lower than the typical GRBs 
studied at high redshifts and it is thus not typical of the 
class. Subsequent modelling of a variety of SNe associ- 
ated with GRBs has revealed a wide distribution in the 
kinetic energy (Ek)^ ejecta mass (M^^) and ^^Ni mas s 
(Mse) parameter space (e.g. 'Nomoto et al."2004", "20061). 
For example, in 2003 SN2003dh/GRB 030329 becam e the 



first bright G RB to have its SN studied in detail ( jHjorth 
et al. | |2003[ [Stanek et al. , 2003) . At redshift z=0.169 



( [Greiner et al. pl)03[ ), it was the first nearby GRB asso- 



ciated with a SN to have 7-ray and afterglow properties 
similar to t hose o bserved for cosmolo gical GRBs ( |Price | 
et al.||2003 ). Also, Berger et al. (2011| recently observec 
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the cosmological GRB/SN 2009nz, and reported broad 
line features with lower velocities than 1998bw. How- 
ever, as explosion parameters are often estimated using 
the simple Arnett relation (Arnett 1982), more detailed 
modelling is warranted (see discussion). 

The collapsar scenario for LGRB-SN predicts that in 
many cases we will be able to observe the SN explosion 
but will not see the GRB because of our orientation. One 
way to test the model is to hunt for the predicted popu- 
lation of asymmetric SN explosions in which our line-of- 
sight is off the jet axis. A promising approach is radio 
surveys looking for emission associated with the inter- 
action of relativistic outflow with a previously emitted 
wind or the c ircumstellar material. B ased on a sample of 
68 Ib/c SNe, iSoderbe rg et al.| ( |2QQ6 ') argued that < 10% 
of type Ib/c SNe are associated with off-axis jets. How- 
ever, radio observations have also shown that outflows 
powered by a central engine may be pr esent in some ob 



jects for which no GRB was detected ( Soderberg et al 



201^. 

Utner phenomena that could be explained as an 
off-axis GRB include X-ray flashes (jYamazaki et al. 



|2QQ4 ), which were also found to sometimes be associ- 



ated with SNe (e.g. XRFlQQ316D/S N2QlQbh; Olivares 
E. et al.||2QT2| iBufano et al.||2QT2| [Cano et al.||i^0111 



XRF Q6Q2 1875N2QQ6aj; |Campana et al^2QQ6 " Plan et al. 



2QQ6[ ). Some evidence lor jet-like SN explosions has 
also been uncovered via late-time optical observations 
of stripped-envelope core collapse SNe. Once the ejecta 
is optically thin the observed blue- and red-shifted lines 
follow the shape of the ejecta. In this stage, double- 
peaked l ine profiles arise wh en a bipolar SN is observed 
off-axis (Maeda et al. 2002), and some cases were de- 



tecte d and analysed by, lorexa mple, Maeda et al. (2008) 
and Taubenberger et al. (2009). 

The aforementioned studies of radio emission and 
late-time optical spectra are powerful tools to search 
for evidence of jets having been involved in SN ex- 
plosions. However, signatures are also expected in 
ultraviolet /optical/near-infrared emission around max- 
imum light of the SN explosion. Observations of this 
phase are typically more complex to interpret than e.g. 
late-phase nebular spectra but are extremely important 
since they are the most routinely accessible data. To 
interpret such data requires suitable explosion models 
and synthetic spectra and light curves. The quest to 
understand the GRB-SN connection has led to the de- 
velopment of numerous hydro dyn amical st udies of jet- 
driven models. Proga et al. ( 20Q3[ ) and later Sawai et al. 
(|2005|) showed via simulations that magnet ohydrody- 



namical processes alone can launch and sustain polar out- 
flows from the accretion disk around a black hole formed 
by the collapse of the core of a massive star. The success 
of such simulations of jet launching have ignited a study 
of theoretical modelling of asymmetric SNe. In particu- 
lar, hydro dynamical simulations of asymmetric explosion 
models have been conducted (Maeda et al. 2002 ; Maeda 
|2006l 'Maeda et al.'"2006l Tanaka et al.j[2007|), and com- 
pared with observations, mostly looking lor a promising 
model to explain the extraordinary SN 1998bw. 

In this paper our goal is to carry out self-consistent, 
multi-dimensional, time-dependent radiative transfer 
simulations to predict the angle-dependent spectra and 
light curves for asymmetric Ic SNe. In previous studies 



(e.g. Tanaka et al. 2007), light curves were calculated 
with a very simplitied treatment of opacity while spec- 
tra were computed separately using a time-independent 
approach. Our new self-consistent calculations will help 
quantify the observable signatures of a jet-driven explo- 
sion and, in particular, allow us to predict the differ- 
ences in optical colors and spectral features that are ex- 
pected for a jet-driven explosion viewed off- axis. As in- 
put, we work with the asymmetric explosion models of 
(Maeda et al.|2002[ , hereafter M02). We also investigate 
the inlLuence ol progenitor metallicity on the optical ob- 
servables since several studies have suggested a possible 
link between the metallicity of the SN envir onment and 
whether it will pro duce a GRB ([Sol lerman et al. 2005 
'M_Qdjaz et al.|2008[ [20061 |Levesqueetalj2010V We begin 



(Q by describing the explosion models and principles of 
our radiative transfer code. In ^we outline our results 
and in ^we discuss our conclusions. 

2. METHODS 

2.1. Explosion Models 

For this study, we adopt the jet-driven explosion mod- 
els of M02. These were created by exp loding the core 
of a s tar that had initial mass of AOMq (Iwamoto et al. 
1998). A cut in mass coordinate (at 2A Mq) was used 
to divide the core model into an inner region, which is 
assumed to collapse to a black hole, and an outer re- 
gion that will form the SN ejecta. 2D hydrodynamical 
and nucleosynthesis calculations were performed to sim- 
ulate the explosion of the outer layers. The explosions 
were driven by inserting an excess of kinetic energy in 
the region below the mass cut. This gives rise to a shock 
which propagates through the model, unbinding the ma- 
terial and triggering nuclear burning. All models have 
a total ejecta mass of ~ 10 M©, although small devia- 
tions arise owing to the treatment of boundary conditions 
in the hydrodynamic simulations. Different asphericities 
were induced in the models by varying the injection of 
the kinetic energy. Specifically, direction-dependent ini- 
tial velocities were imposed such that the velocity along 
the polar (i.e. assumed jet) axis (v^ = Oiz) was larger 
than that along the equatorial direction {vr = Pr). The 
degree of asphericity is conveniently parametrized by the 
ratio a/f3. 

The M02 models explored a range of values for the fi- 
nal kinetic energy and degree of asphericity. Since this 
study focuses on the effect of asphericity on the light 
curves and spectra, we will consider only a subset of their 
models that have a fixed final kinetic energy of ^51 = 20 
(^51 = Eiso/lO^^eigs): models A, C and F of M02. Of 
these, model A is the most aspherical, C has intermediate 
sphericity and model F is spherically symmetric. Impor- 
tant parameters of these models, and their ^^Ni yields, 
are given in Table [ij 

Observations of SNe associated with LGRBs do not 
reveal any H or He. Modelling of SN1998bw was suc- 
cessful when using as a progenitor a core stripped down 
to its O layer by the time of explosion (Maeda et al.| 
2002). Therefore, for our radiative transfer calculations 
the composition of all unburnt material has been set to 
that of the O layer with a small fraction (~1%) of C (i.e. 
we do not include any He-rich material in the ejecta). 
To explore the effect of progenitor metallicity, we have 
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TABLE 1 
Explosion Models 



Model 


Mej{MQ) 








Al (AO.l) 


10.10 


0.23 


20 


16 


CI (CO.l) 


10.05 


0.28 


20 


8 


Fl (FO.l) 


10.08 


0.38 


20 


1 



^ Asphericity parameter — see text for details 

included all abundances of elements heavier than Na (up 
to atomic number 30) in the composition of the unburnt 
material. We adopt ed the solar elemental abundances of 
[Anders fc Grevesse] ( 1989) and have computed synthetic 
observables tor models with metallicities (relative to so- 
lar) of Z = 0.1 Zq and Z = IZq. Note that the solar 
element abundances have only been added to the input 
conditions for our radiative transfer simulations - i.e. we 
explore the effect of metallicity in the unburnt material 
on the synthetic spectra and light curves but do not in- 
vestigate how progenitor metallicity influences either the 
hydrodynamics or nucleosynthesis of the burnt material. 
In addition, the metallicity of the primordial material 
of a s tar will influence its evolution (e.g. Meynet et al. 



2009). This will effect the abundances and stellar struc 
ture at the time of collapse, which will alter the expected 
SN spectra and light curve. We plan to investigate the 
expected variations in a coming paper (Rapoport et al. 
in prep) in which we consider a wider range of progenitor 
models. Throughout this paper, we will denote our cal- 
culations adopting solar metallicity for the M02 models 
A (C, F) as Al (CI, Fl) and those for sub-solar metal- 
licity {Z = O.lZo) as AO.l (CO.l, FO.l). 

Fig. [l] illustrates the spatial distribution of mass den- 
sity and the abundances for key elements that will feature 
in our discussion. In the aspherical models (Al, CI), the 
explosion dynamics gives rise to a relatively dense core 
region that is extended in the equatorial direction. In 
contrast, the polar region (i.e. the direction in which the 
high velocities were initially injected) has relatively low 
density. ^^Ni is produced by explosive Si burning and, in 
all aspherical models, its mass fraction is highest in the 
low density regions. This gives a near-conical shape to 
the ^^Ni distri bution for the most aspherical model (see 
also Fig. 2 in Maeda et al. 2006). When the degree of 
asphericity is reduced (model CI), the ^^Ni distribution 
no longer flares out from the pole and gradually looks 
more like the spherically symmetric model (Fl). 

Si (and other intermediate-mass elements) is produced 
by explosive O burning (e.g. ^Nomoto et al.,,2006) . Its 
abundance distribution is similar in shape to ^'^Ni, being 
concentrated around the outer edge of the ^^Ni-rich re- 
gion. Si is also present in small quantities in the outer 
layers since we have assumed non-zero progenitor metal- 
licities in all models. The Na in our models is domi- 
nated by the contribution in the unburnt material. Its 
abundance distribution is therefore comparatively uni- 
form and its concentration highest in the outer layers. 



2.2. Radiative Transfer 
Our radiative transfer simulations were carried out us- 
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Fig. 1. — Left column: 
right column: Model Fl. 



Model Al, middle column: Model CI, 
From top to bottom: Total mass den- 



sity, ^^Ni, Si and Na mass fractions. While model F is spherically 
symmetric the 'line' shape seen for the Si and ^^Ni mass fractions 
is due to numeric effect as hydrodynamic simulations tend to give 
artefacts near the imposed axes of symmetry. The bulging along 
the equatorial is another simulation feature due to numeric effects 
and the colliding at the equator of the bow shocks which are driven 
by the jet. 

Monte Carlo (MC), multi-dimensional, time-dependent, 
radiative transfer code based on the methods described 
2005). The code computes syn- 



by Lucy (2002 2003 



ing ARTIS (Sim^,2007| [Kromer fc Sim^2009j , which is a 



thetic spectra tor explosion models assuming homol- 
ogous expansion of the ejecta and the Sobolev approxi- 
mation for line opacity. MC radiative transfer methods 
has been applied to a variety of SN explosion mod els and 
tested by co mparison to other calculations in e.g. |Kasen 
et al l (|2QQ6t ) (using the SEDONA code) and [KromerT 
5im| l [20ngi r(for ARTIS). Here, for the first time, we ap- 
ply 177775 to models for SNe Ic. We note that, for 
SNe Ic explosions with a compact progenitor, homolo- 
gous expansion is expected to be reached within the first 
~10 seconds and is therefore a good approximation at 
the times of ~days that will be studied here.). 

As input, the code requires specification of the veloc- 
ity, density and composition for each grid cell (a Carte- 
sian grid is used). The total energy emitted from the 
decay chain ^^Ni^^^Co^^^Fe is determined from the 
model and is divided into a set of MC quanta. At the 
start of the radiative transfer simulation, these quanta 
are placed in the model according to the ^^Ni densities 
in the form of radioactive pellets (Lucy 2005). During 
the simulation these convert into photon packets in ac- 
cordance with the radioactive decay times. The photon 
packets then propagate through the model of the ejecta 
and can interact with the medium via Compton scatter- 
ing, pair production, bound-free and bound-bound ab- 
sorption. For bound-bound transitions we use a line list 
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TABLE 2 

V - BAND MAXIMUM LIGHT TIMES IN 
DAYS AFTER EXPLOSION 



V 



Model 25° 



45° 60° 85° 



Al 


22 


24 


24 


24 


AO.l 


20 


22 


22 


22 


CI 


24 


24 


24 


24 


CO.l 


23 


24 


24 


24 


Fl 




29 






FO.l 




28 







^ Values are quoted for four observer 
inclination angles (measured relative 
to the polar axis). 

^ One number for all columns implies 
no angle dependency. 

extracted fro m the dat a of "Kurucz fc Bell (1995") CD23, 
as describe in'Kromer & Sim (2009). The time, direction 
and frequency of escaping energy-quanta are recorded to 
allow angle-dependent spectr a and li^ht curves to be cre- 



ated. We refer the reader to Kromer & Sim (|2009) for 
full details on the MC techniques employed by ART IS. 

Spectra were generated for models Al, AO.l, CI, CO.l, 
Fl and FO.l. The simulations were started at three days 
after explosion and were run until 50 days after explo- 
sion with 60 time steps, logarithmically spaced. While 
a study of the models at earlier start time is appealing, 
the high opacities immediately after explosion results in a 
heavy computation requirement. In addition, this paper 
focuses on SNe associated with real GRBs, where early 
time SN properties are practically undetectable under 
the luminous GRB afterglow - consequently predictions 
for very early times are not easily testable. A grid of 
68x68x68 cells was used, giving a resolution of ~ 1, 180 
km/s. As the explosion models are two dimensional with 
symmetry along the jet direction, the synthetic observ- 
ables depend only on the angle {0) between the observer 
line-of-sight and the polar axis of the model. For the as- 
pherical models, we extracted synthetic observables for 
5 values of by dividing the emergent MC quanta into 5 
bins of equal solid angle across the range < 6> < 90 deg. 

3. RESULTS 
3.1. Light curves 
The main advantage of our simulations com pared to 



et 



previous stu dies of the M02 explosion models (Tanaka 
al. 12007 hereafter TOT) is a self-consistent 
)ena( 



time- 

dependent treatment of the spectral evolution. This 
means we can extract light curves that are fully con- 
sistent with both the spatial distribution of ^^Ni in the 
model and the detailed opacities needed to compute re- 
alistic s pectra. Fig. [ 2 | shows ligh t curves for the standard 
Besseh ^Besseh fc Miirphy|2012[) U, B, V, R and / filterf] 



for Model Al lor lour ditterent observer orientations. We 
also show the bolometric UVOIR light curves, defined 
by a top-hat filter extending from 3,000A to 9,000A. 

In 5, V ^ R and /, the light curves for all observer 
orientations reach peak around 20 - 25 days after ex- 
plosion, (~7 days later than typical for SNe associated 

^ For convenience, the filter functions are shown in the upper 
panel of Fig. [s] 



with GRBs, see discussion) and then decline monoton- 
ically. In these bands, the time of peak depends only 
weakly on the observer orientation but tends to be earlier 
for inclinations close to the pole (small ^), particularly 
in the bluer bands. In U band, the infiuence of orien- 
tation is much stronger - the U light curve rises much 
more quickly and peaks around ten days earlier for a po- 
lar compared to equatorial inclination. The widely-used 
light curve decline rate parameter Amis (the increase in 
magnitude from peak to 15 days after peak), is similar 
for the asymmetric models and vary slightly with angles 
(see Table [3|. 

In general, the pre-maximum light curves are brighter 
for polar inclinations but around peak the equatorial 
light curves become brighter. Since the angle depen- 
dence is strongest in the bluest bands, there are also sig- 
nificant differences in the color evolution for polar and 
equatorial lines of sight. This is illustrated in Fig. [3j 
which shows the B — V color evolution for Model Al. 
If viewed close to the polar direction, the color becomes 
dramatically redder during the evolution towards max- 
imum light. In contrast, the equatorial color evolution 
is much weaker, showing increasing blue- ward evolution 
starting a few days prior to maximum light. Similar 
orientation-dependent trends are predicted in our other 
asymmetric models. The scale is smaller in the more 
spherical models (in Model CI the difference between 
^ — V for polar and equatorial viewing angle at max- 
imum light is 0.1 mag.) and larger in corresponding 
models with reduced metallicity (e.g. Model AO.l, where 
the difference is 



3.4). 



0.45 mag.; see Section 
The origin of the angle-dependence in our light curves 
and colors can be largely traced to the asymmetric dis- 
tribution of iron group nuclei, particularly ^^Ni and its 
decay products (see Fig. [T]) , and the evolution of the ion- 
ization state of the ejecta, which is shown in Fig. [4j At 
early times, the effective photosphere is located at rela- 
tively high velocities. Since the ^^Ni distribution is most 
extended around the poles, it is there that the photo- 
sphere is first most directly heated, leading to initially 
stronger polar emission. The ionisation state is gener- 
ally high at early times (e.g. in the ^^Ni-rich regions, Co 
is dominated by Co ill as seen in Fig. |4][^ ), such that 
this emission is fairly blue. As the ejecta expands, the 
opacity of the outer ejecta drops and the photosphere 
recedes deeper into the ejecta. Once radiation starts to 
escape from the inner ejecta, the larger projected area 
of the ^^Ni-rich region for equatorial inclinations gener- 
ally favours brighter peak magnitudes for these orienta- 
tions. At the same time, the line blanketing in the blue 
is stronger along the poles where metals synthesised in 
the explosion are present in the region around the pho- 
tosphere. This leads to relatively red colors. This iron- 
group material is mostly singly- and doubly- ionised (up- 
per panels of Fig. B but it gradually recombines, causing 
the polar line blanketing to become even stronger with 
time in the rise to maximum light. 



^ In Fig. ^ we chose to show Co as an indicator for the typical 
ionization conditions for iron group elements synthezised in the ex- 
plosion. For the epochs shown in the plot, ^^Co is the most abun- 
dant iron-group isotope although we note that other iron group 
elements (including Ti, Cr and Ni) are still very important for line 
blocking in the blue. 
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TABLE 3 

Amis - THE INCREASE IN MAGNITUDE FROM PEAK TO 15 DAYS AFTER PEAK IN [/, B, V, R AND / BANDS FOR DIFFERENT OBSERVER INCLINATION 

ANGLES 



Model 




U 








B 






V 






R 






/ 




25° 1 


45° 


60° 


85° 


25° 


45° 60° 


85° 


25° 


45° 60° 


85° 


25° 


45° 60° 


85° 


25° 


45° 60° 


85° 


Al 


0.4 


1.0 


1.0 


0.8 




0.72 






0.9 






0.9 




0.7 


0.6 0.7 


0.6 


AO.l 


0.6 


1.0 


1.1 


1.0 




1.1 






1.0 






1.0 






0.6 




CI 


0.9 


0.7 


0.9 


0.9 




0.7 




0.9 


0.8 0.8 


0.8 


0.9 


0.9 0.9 


0.7 


0.7 


0.6 0.6 


0.6 


CO.l 




1.0 






1.1 


1.0 1.1 


1.0 


1.0 


0.9 0.8 


0.9 


0.9 


0.9 0.9 


0.7 


0.6 


0.6 0.6 


0.5 


Fl 




1.0 








0.6 






0.6 






0.6 






0.4 




FO.l 




1.2 








1.0 






0.7 






0.6 






0.4 





^ Values are quoted for four observer inclination angles (measured relative to the polar axis). 
2 One number for all columns implies no angle dependency. 



-19 



Model Al 




10 15 20 25 30 
Time after explosion (days) 

Fig. 2. — Light curves for model Al for the different viewing 
angles. ^ = 0° is the jet direction and = 90° is the equator. See 
top panel of Fig. [sjfor the filter curves. 



3.2. Si II velocity and correlation with color 

The Si II 6355A transition is an unmistakable absorp- 
tion feature that is commonly seen in the red part of SN 
spectra. The asphericity of the Si-rich region (see Fig. 
[1]) gives rise to a significant angle-dependence of the Si ll 
6355A absorption line velocity. This was previously re- 
ported by T07 and is confirmed in our study (see Fig. [sj 
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Fig. 3. — B-V color evolution for model Al for the different view- 
ing angles. ^ = 0° is the jet direction and = 90° is the equator. 
See top panel of Fig. [sjfor the filter curves. 




10 20 30 10 20 30 10 20 30 10 20 30 
r velocity (10^ km/s) 

Fig. 4. — Ion fractions for model Al. From left to right the ions 
increase from I to IV, from top to bottom the ions are Co at 12 
days (~10 days before V-Band maximum lights), and Co, Si and 
Na at y-Band maximum light (see table |2] for details). The blue 
curves are contours of the ^^Ni mass fraction and correspond to 
0.05, 0.07 and 0.09. 
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Fig. 5. — Synthetic spectra at V-band maximum (see table [2] for 
details) light for models Al (most asymmetric), CI (intermediate 
symmetry) and Fl (spherically symmetric, upper, middle and lower 
panels, respectively) for different viewing angles. The top panel 
shows the filter functions used to generate the light curves in Fig. 
[2] ^ = 0° is the jet direction and = 90° is the equator. 



which shows the optical spectra computed for our models 
at F-band maximum light). To illustrate the origin of the 
angle-dependence of the Si line, Fig. [6] shows the Sobolev 
optical depth of the 6355 A line at 15, 25 and 35 days for 
models Al, CI and Fl and indicates the location of the 
region of last interaction for escaping i?-band photons. 
At all epochs, the region in which the line is optically 
thick extends to significantly higher velocity along the 
polar than equatorial direction (in models Al and CI). 
Therefore the absorption extends to higher velocity in 
the spectrum when observed pole-on. Notice, however, 
that in most cases the region of last interaction for es- 
caping i?-band photons does extend outside the region in 
which the Si line is optically thick. Consequently, the Si 
line trough has usually been partially refilled such that 
the line core does not appear to be saturated. 

The same trend appears across the spectrum and can 
clearly be seen in the Ca ll triplet, especially on the blue 
side of the emission part. Higher velocities towards the 
equator are also seen in the blue (e.g.~4400A), although 
the overlapping of many Fe, Ti and Co lines makes quan- 
titative analysis more complicated in this region. 

As discussed in Section [^T] the SN colors at maximum 
light also show a systematic dependence on the observer 
orientation. Consequently, our simulations predict that 
the Si velocity and the color should be correlated for SNe 
described by the M02 models. This predicted correlation 
is shown in Fig. [7| For all our aspherical models, we 
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Fig. 6. — Sobolev optical depth (color coded) of the Si 11 6,355A 
line for models Al, CI and Fl (from top to bottom) at times 15, 25 
and 35 days after explosion (from left to right). The black contours 
indicate the area of last interaction of escaping R band photons at 
that time (effective photosphere). The contours correspond to 0.2, 
0.5 and 0.8 of the emissivity normalised to peak value. 

find that decreasing the inclination angle causes the peak 
B — V color to become systematically redder while the 
Si II velocity becomes simultaneously higher. The range 
of variation depends on the degree of asphericity. Clearly, 
the spherically symmetric models show no variation with 
observer orientation while in Model Al we find the B-V 
color to differ by up to ~0.16 magnitudes and the range 
of Si II velocities (A^si n) to be -6,300 km s-^ Model CI 
shows a similar correlation to that in Model Al, except 
that the range of variation is smaller since the degree of 
asphericity is lower. 

We find that this correlation between line velocity and 
color is present in the simulations from —18 days after ex- 
plosion for as long as the Si line remains in the spectrum. 
However, it changes slightly at different epochs. For ex- 
ample, in Model Al, the absolute Si velocities decrease 
after maximum light and the amplitude of the color vari- 
ation with viewing angle becomes slightly smaller (see 
Fig.[3|. Nevertheless, the correlation persists and, for 18 
days after explosion and all later epochs, is always in the 
same sense (redder colors are associated with higher line 
velocity) . 

3.3. Na I absorption 

One of the strongest features in our synthetic spectra 
is the Na I 5890A absorption line (see Fig. |5|. In gen- 
eral, we find this feature to have similar strength and 
angle dependency as reported by TOT. The Na in the 
ejecta is predominantly in the unburned material. Its 
abundance is therefore governed by the Na mass-fraction 
in the pre-explosion stellar model and is uniform in the 
outer ejecta (Fig. T]). Due to the high optical depth of 
this ground state transition, only a small Na mass frac- 
tion is required to make the line strong, if the ionization 
state is sufficiently low to favour Na I. The angle de- 
pendence of the Na line is therefore not a consequence 
of the spatial distribution of the Na abundance, but of 
variations in the ionisation state (see Fig. [8|. Particu- 
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Fig. 7. — B — V color vs. Si ii velocity for all the models at 
y-Band maximum light (see table [2] for details). For the aspheri- 
cal models (Al, AO.l, CI, CO.l) the symbols represent data points at 
25,45,60,72 and 85 degrees with decreasing Si velocity (from right 
to left). The horizontal error bars indicate the uncertainties asso- 
ciated with fitting the line velocity from our MC spectra. The lines 
connect the data points to guide the eye. 
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Fig. 8. — Na l (left) and Na 11 (right) density distribution for 
model Al at V-Band maximum light (see table [2] for details). In 
both panels, the distributions are normalized to a maximum value 
of 1.0. The neutral Na is concentrated along the ~ 90° direc- 
tion explaining the deep absorption observed from an equatorial 
orientation. 



larly for early times, it is only in the relatively dense 
equatorial regions that the Na I population is sufficiently 
high that the line is very strong (see Fig.[9|. However, at 
later times sufficient recombination occurs that the line 
becomes significant for all orientations. 

As noted by T07, the GRB related SN 1998bw did not 
show strong Na I absorption, which is consistent with 
the models for a pole-on viewing angle around maximum 
light. However, the Na I line remains a challenge for our 
models since the synthetic observables for equatorial ori- 
entations are expected to correspond to cases in which 
the jet axis is not close to our line-of-sight. Such off-axis 
explosions should produce a sub-population of SNe Ic. 
However, strong Na I is not typically observed for broad- 
lined core collapse SNe. This may simply be a failing of 
the particular explosion models and/or ionization treat- 
ment adopted here. As discussed above, the formation of 
this feature is very sensitive to the degree of ionization in 
the outer ejecta and it could be suppressed if the ioniza- 
tion remained slightly higher. Increased ionisation could 
be achieved in models with somewhat higher ^^Ni masses 
(as woul d be r equired to account for the brightness of e.g. 
1998bw (|M^da 2006)). 
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Fig. 9. — Sobolev optical depth (color coded) of the Na i 
5,890A line for models Al, CI and Fl (from top to bottom) at 
times 15, 25 and 35 days after explosion (from left to right). The 
black contours indicate the area of last interaction of escaping V 
band photons at that time (efective photosphere). The contours 
correspond to 0.2, 0.5 and 0.8 of the emissivity normalised to peak 
value. 
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Fig. 10. — Synthetic spectra for models Al and AO.l at V-band 
maximum light for polar (r-^ 25°) and equatorial (~ 85°) viewing 
angles. The top panel shows the filter functions used to generate 
the light curves in Figs. [2] and [n] 



3.4. Metallicity Dependence 

Studies reveal a possible link between the metallicity 
of the SN environment and whether it will produce a 




by the small sample of spectroscopically 
firmed GRB- SNe, only one of which has had sufficiently 
complete observations to be modelled in detail (2003dh, 
Mazzah et al.""20031 'Matheson et al."'20031 'Stanek et al. 



2003; Simon et al. 2004; D eng et al. 2005) . Sollerman 
et al ( 2005[ ) observed three galaxies which were identi- 
tied to host a SN with an associated GRB. Their results 
favour low metallicity, sub-luminous galaxies in a phase 
of active star formation, whi ch agrees with ot her stud- 
ies of higher redshift GRBs (iLe Floc'h et al.||2003). It 
is therefore important to consider how our synthetic ob- 
servables are affected by a reduced progenitor metallicity, 
to confirm the robustness of our findings and identify po- 
tential observable signatures of progenitor metallicity in 
the SN itself. 
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Fig. 11. — Light curves for model AO.l for the different viewing 
angles. ^ = 0° is the jet direction and 6 = 90° is the equator. The 
shaded gray area spans the values for all angles for model Al to 
allow comparison. See top panel of Fig. [sjfor the filter curves. 



Fig. 10 shows spectra for our sub-solar metallicity 



model AO.l and solar metallicity model Al for the same 
epoch as shown for model Al in Fig. [5j The biggest dif- 
ference is increased flux in the blue part of the spectrum 
for all orientations, a consequence of less line blocking 
from primordial heavy elements in the lower metallicity 
model. This persists for all epochs and results in light 
curves that are systematically brighter up to ~35 days 
past explosion and reach peak earlier, in B and U for 
model AO.l (Fig. 11) compared to model Al. We draw 
similar conclusions when comparing the spherically sym- 
metric models. 

Since the metallicity affects the redder bands much less 
significantly, the increased blue flux from the sub-solar 
metallicity models influences the colors. The sub-solar 
models are always bluer and also show a wider varia- 
tion of color with orientation. This makes the Si ll-color 
correlation even stronger (and, consequently, more de- 
tectable), spanning a range of A{B — V) ~0.5 around 
peak (see results for model Al in Fig. I7|. The trend 
in the correlation with time is different to that found in 
model Al with a plateau around peak and a beginning 



of decrease ~35 past explosion. 

4. DISCUSSION/CONCLUSION 

We performed multi-dimensional radiative transfer cal- 
culations to compute light curves and spectra for 2D 
models of SNe driven by bi-polar outflows. Compared 
to previous studies, the main advantage of our work is 
the use of self-consistent opacities in a time-dependent 
calculation of the spectral evolution. 

As expected for strongly asymmetric explosion mod- 
els, we find that many observable properties depend on 
the orientation of the observer. During the rise phase, 
the synthetic light curves are brighter and the colors are 
bluer for an observer inclination close to the polar axis. 
However, at peak, the pole-on light curves are slightly 
fainter and redder than those seen from an equatorial 
orientation. This is a consequence of the concentration 
of ^^Ni along the polar direction in the explosion mod- 
els. The time of maximum light also depends slightly on 
inclination and is reached a few days earlier for pole-on 
inclinations, in agreement with the findings by [Maeda 
( [20061 ). 

T he color evolution of our sub-solar metallicity aspher- 
ical model AO.l in the pole-on case shows similar trends 
to 1998bw: B — V increases with time, although less 
quickly in our model than for 1998bw. We note, how- 
ever, that 1998bw is systematically bluer at all times by 
0.2-0.4 magnitudes. In our models, we find low metallic- 
ity makes the peak B — V color bluer and increases the 
steepness of the B — V temporal color evolution around 
peak, suggesting that even lower metallicity might be 
worthy of investigation for 1998bw. Comparing our off- 
axis color evolution to SN 2009bb, a BL-Ic SN without 
a GRB, we find similar results (i.e. with a slower color 
ev olution and redder col or in our model than those found 
by|Pignata et al.|(|201l| for the SN). 



Our asymmetric models predict that an on-axis orien- 
tation gives rise to a redder B — V color at peak. This 
should be statistically testable if a large sample of poten- 
tial jet-driven SNe Ic (including events observed both on- 
and off-axis) were identified. This could be most clearly 
done with a set of objects for which the ejecta mass and 
energy are similar. However, even for a sample with a 
range of explosion energies and masses a mean trend in 
the ensemble averages might be detectable. 

These orientation-dependent effects on the light curve 
mean that one would overestimate the ejected mass and 
^^Ni (by on the order of tens per cent) if a pplying 
the widely-used Arnett relation (Arnett 1982) to the 
equatorially- viewed light curve, and underestimate these 
parameters when viewing the SN pole-on. 

Spectral features also depend on the viewing angle. As 
already noted by TOT, the Si ll line velocity at maximum 
light varies monotonically (within our MC uncertainties) 
with orientation, being fastest for pole-on inclinations. 
We find the Si ll velocities along the poles to be compa- 
rable with those reported for SNe associated with GRBs 
such as 1998bw and 20031w, while our off- axis velocities 
are comparable to BL-Ic SNe such as 1997ef and 2003bg 
( Corsi et al.||2011 ). Since our light curve properties also 
show a simple trend with viewing angle, this means we 
predict a Si ll velocity - peak color {B — V) correla- 
tion. We found that this correlation is robust: it is pre- 
dicted for solar and sub-solar models during a significant 
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range of observable epochs and for a range in aspherici- 
ties. Such a correlation should be statistically testable if 
a large sample of potential jet-driven SNe Ic (including 
events observed both on- and off-axis) were identified. 

We examined the role of progenitor composition on the 
light curves and spectra by considering models appropri- 
ate for both solar and one-tenth solar metallicity. We 
showed that the metallicity strongly infiuences the blue 
part of the spectrum (the B — V color is roughly 0.2 - 
0.8 mag bluer for Z = O.IZq compared to Z = IZq). 
Therefore, if GRB progenitors typically have low metal- 
licity, we expect them to have measurably bluer colors 
than a typical stripped-envelope core collapse SN with 
a more metal-rich progenitor. However, the metallicity 
does not qualitatively affect the trends of line-velocity 
and light curve color with orientation (indeed, we find 
that our velocity-color relationship is even stronger for 
Z = O.IZq). 

There are some shortcomings of our synthetic light 
curves and spectra in comparison to observations. Typ- 
ically, our light curves reach peak at ~23 days, which 
is 7 days lat er than typ ical for SNe associated with 
GRBs (Woosley & Bloom 2006) which could imply that 
the density distribution in the current models is some- 
what inappropriate. Also, at most epochs in our models, 
we predict strong Na I absorption, which is not typi- 
cal of SNe Ic. This may be attributed to the ionisation 
state being too low in our current models. Further study 
and exploration of different explosion model parameters 
is warranted to understand whether these issues can be 
resolved. 

In this work we used the the well-studied M02 mod- 
els as a test case to compare our methods with previ- 
ous results, and allow us to understand the limitations 



and strengths of our approach before applying them to 
a wider range of models. We have concentrated on only 
two of the parameters that are relevant to asymmetric 
Type Ic explosions: the degree of asphericity and the 
composition (i.e. metallicity). However, to fully explore 
this class of models and make quantitative comparisons 
to observed explosions, we will extend this work to span 
a range of masses and explosion energies. Having such a 
grid of synthetic spectra for models with different explo- 
sion parameters will help to quickly identify parameter 
ranges for newly observed GRB- SNe, on which more de- 
tailed modelling can be based. 
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